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Nonresonance and global existence of 
prestressed nonlinear elastic waves 

By Thomas C. Sideris* 
1. Introduction 

This article considers the existence of global classical solutions to the 
Cauchy problem in nonlinear elastodynamics. The unbounded elastic medium 
is assumed to be homogeneous, isotropic, and hyperelastic. As in the theory of 
3D nonlinear wave equations in three space dimensions, global existence hinges 
on two basic assumptions. First, the initial deformation must be a small dis- 
placement from equilibrium, in this case a prestressed homogeneous dilation of 
the reference configuration, and equally important, the nonlinear terms must 
obey a type of nonresonance or null condition. The omission of either of these 
assumptions can lead to the breakdown of solutions in finite time. In particu- 
lar, nonresonance complements the genuine nonlinearity condition of F. John, 
under which arbitrarily small spherically symmetric displacements develop sin- 
gularities (although one expects this to carry over to the nonsymmetric case, as 
well), [4]. John also showed that small solutions exist almost globally [5] (see 
also [10]). Formation of singularities for large displacements was illustrated by 
Tahvildar-Zadeh [16]. 

The nonresonance condition introduced here represents a substantial im- 
provement over our previous work on this topic [13]. To explain the difference 
roughly, our earlier version of the null condition forced the cancellation of 
all nonlinear wave interactions to first order along the characteristic cones. 
Here, only the cancellation of nonlinear wave interactions among individual 
wave families is required. The difficulty in realizing this weaker version is 
that the decomposition of elastic waves into their longitudinal and transverse 
components involves the nonlocal Helmholtz projection, which is ill-suited to 
nonlinear analysis. However, our decay estimates make clear that only the 
leading contribution of the resonant interactions along the characteristic cones 
is potentially dangerous, and this permits the usage of approximate local de- 
compositions. 
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The limitations of the earlier version of the null condition were exposed in 
the work of Tahvildar-Zadeh who first considered small perturbations from an 
arbitrary prestressed homogeneous dilation [16]. Using an expansion based on 
small shear strains, he observed that the null condition in [13] placed physically 
unrealistic restrictions on the growth of the stored energy function for extreme 
values of the dilational strain. Our revised null condition avoids this defect. 
In fact, an expansion involving small distortional strain (which to first order 
is equivalent with small shear strain) reveals that the dominant behavior of 
the stored energy function for arbitrary dilational strain is determined by the 
propagation speeds of the medium. The null condition enters as a higher order 
correction. 

Outside of the inevitable energy estimates, the techniques of the existence 
proof represent an evolution of the ideas initiated in [10] and expanded in [13]. 
The central hurdle to be overcome is the fact that the equations of motion for 
elasticity are not Lorentz invariant. For the nonlinear wave equation, a well- 
oiled machine based on Lorentz invariance, the so-called generalized energy 
method, has been developed by Klainerman for the construction of solutions [8]. 
Its attractive feature is the avoidance of direct estimation of the fundamental 
solution of the wave equation through the use of weighted Sobolev inequalities 
involving the generators of the Lorentz group. The generators are effective 
because they commute with the wave operator. In the case of elasticity, it 
is still possible to get weighted estimates with respect to a smaller number 
of generators, but the resulting estimates are weaker. In particular, the best 
decay is available only for second (or higher) derivatives of the solution. In our 
earlier work [10], [13], this was dealt with by smoothing and the consequent 
introduction of a somewhat unnatural and cumbersome nonlocal energy. The 
proof presented here represents a refinement also in that the use of the nonlocal 
energy is now completely avoided. Given the significant differences from our 
earlier work, we have endeavored to make a self-contained presentation. 

The results presented herein were announced in [14]. R. Agemi has re- 
cently completed a manuscript dealing with the existence of solutions near the 
unstressed reference configuration under the same null condition used here [1]. 
His proof relies on direct estimation of the fundamental solution. The approach 
taken here can be adapted to the case of a system of coupled quadratically non- 
linear wave equations in 3D with multiple propagation speeds [15]. 

The remainder of this introduction will be devoted to the description of 
the basic notation and the formulation of the equations of motion, leading as 
speedily as possible to a statement of the global existence theorem. In the 
following section, we then explore the origins of our assumptions in an effort 
to convince the reader of their transparency. The rest of the paper is devoted 
to the proof of global existence. The key components: commutation, null form 



NONLINEAR ELASTIC WAVES 



851 



estimates, weighted L°° estimates, and weighted 1? estimates, are assembled 
along the way in a series of four propositions. 



1.1. Notation. Partial derivatives will be written as 

— and Oi = — —, 

8t ox 1 

We will also abbreviate 



do = dt 



1,2,3. 



d = (80,81,82,63), and V = (61,82,63). 

The so-called angular momentum operators are the vector fields 

= (0 1 ,0 2 ,fi 3 ) = iAV, 

A being the usual vector cross product. Crucially, the spatial partial derivatives 
can be decomposed into radial and angular components 



1.1] 



V 



-8 r 



A f2, where 



and d r 



V. 
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A central role is played by the generators of simultaneous rotations which 
will be seen in Section 3.1 to arise naturally in connection with the symmetries 
of isotropic materials. They are given by the following vectorial modification 
of the angular momentum operators 



(1.2a) 
with 
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Equally important will be the scaling operator 

S = tdt + rd r , 

but we also shall show in Section 3.1 why it is more precise to use the generator 
of changes of scale 



(1.3) 



S = 5-1. 



The eight vector fields will be written as T = (To, . . . , Tr) = (8, 0, S). The 
commutator of any two T's is either or is in the span of T, and in particular, 
the commutator of V and T is in the span of V. Schematically, we write 

(1.4) [v,r] = v. 

By T a , a = (ai, . . . , a K ), we denote an ordered product of k = \a\ vector fields 

r -r 

1 ai 1 a K - u 



1 Notc that this differs from the standard multi-index notation. 
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The energy associated to the linearized operator to be defined in Section 
1.2 is 

Ei(u{t)) = \ f \\d t u(t)\ 2 + c 2 \Vu(t)\ 2 + (c 2 - c|)(V • u(t)) 2 } dx, 
and higher order energies are defined through 

E K (u(t))= J2 Mr a <t))- 

\a\<K-l 

In order to describe the solution space we also introduce the time-independent 
analog of T. Set 

A=(Ai,...,A 7 ) = (V,fi,r$.-l). 
Then the A's have the same commutation properties as the T's. Define 
HI = {/ € L 2 (R 3 ) 3 : A a f £ L 2 (R 3 ) 3 , \a\ < «}, 

with the norm 



\ = E n Aa /ii| 2( 



\a\<K 

The solution will be constructed in the space Hfi(T) obtained by closing 
the set C°°([0,T);C^(M 3 ) 3 ) in the norm sup El' 2 {u{t)). Thus, 

0<t<T 

H£(T) c\u(t,x) : du(t,-) G f]' C'([0, T); H^) 

[ 3=0 

By (3.20a), it will follow that H£(T) C C K ~ 2 ([0,T) x R 3 ). 

We define the orthogonal projections onto radial and transverse directions 

by 

(iy rp 1* I If 

1.5a) Piu(x) = — (8> —u(x) = — (—, uyx) 

r r r \r 

and 

(1.5b) P 2 u{x) = [I- Pi]u{x) = — — A (— A u(af 

r \r 

Employing the notation (p) = (1 + I/?] 2 ) 1 / 2 , we will use the following weighted 
L 2 -norm 

2 3 3 

(1.6) W)) = EEE E ||<c Q i-r)P a ^r a u(i)|| L2 . 

a=l/3=0£=l |a|<«-2 

1.2. The equations of motion. Consider a homogeneous elastic material 
filling space. Assume that its density in its undeformed state is unity. The fun- 
damental unknown is the deformation ip : R x R 3 — > R 3 which is an orientation- 
preserving diffeomorphism taking a material point x £ R 3 in the reference con- 
figuration to its position tp(t,x) G R 3 at time t. Consequently, the deformation 
gradient F = Vip must satisfy det F > 0. 
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A hyperelastic material is one for which there exists a stored energy func- 
tion cr(F), representing the potential energy term in the Hamiltonian. The 
Piola-Kirchhoff stress tensor is given by da(F)/dF. A material is frame indif- 
ferent, respectively, isotropic if the conditions 

(1.7) e{QF) = a(F) and a(FQ) = a(F) 

hold for every proper orthogonal matrix Q. It is well-known that (1.7) implies 
that the stored energy function a depends on F only through the principal 
invariants of the strain tensor F T F; see [11, p. 192]. 

Under these assumptions, the equations of motion are! 

d 2 ip l d I da 



li 



dt 2 dx l \ dF\ 



0. 



We shall consider deformations of the form </?(t, x) = Xx + u(t,x), A > 0, in 
which u(t, x) represents a small displacement from a homogeneous dilation, 
so that, in particular, ip(t, •) is a bijection. For small displacements u, the 
standard linear theory outlined in Section 2.2 ensures that the system (1.8) 
is hyperbolic. The long-time behavior of solutions of quasilinear wave equa- 
tions in 3D is determined by the structure of the quadratic portion on the 
nonlinearity [2], [3], [7], [9], [12], and so there is no essential loss of generality 
in considering the truncated version of (1.8) obtained by expansion to second 
order about the equilibrium. This yields the quasilinear system 

d 2 u — Au = N(u, u), 



(1.9a) 
with 



(1.9b) 



and 



(1.9c) 



Lu 

(A u y = Ay m (\)d e d m u^ 

8 2 a 



dFfdFi 



■(A/), 



N(u,vY = B 



B emnW 



ijk 
imn 



(\)d e (dr, 



8 3 a 



OFidFLdFii 



(XI). 



When it is secondary, the dependence of the coefficients on the parameter A will 
be suppressed. It follows from the definitions (1.9b), (1.9c) that the coefficients 
are symmetric with respect to pairs of indices 



(1.10) 



A %3 



A 31 



and B 



ijk 
imn 



B 



jik 
mln 



B 



ikj 
inm' 



2 Repeated indices are always summed, regardless of their position up or down. 
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Thus, the form (1.9c) is symmetric and 

(1.11) N(u,v) = N(v,u). 

1.3. The global existence theorem. The basic assumptions of linear elastic- 
ity, summarized in Section 2.2, ensure the hyperbolicity of the linear operator. 
It will be seen in (2.6b) to have the form 

(H) L = (d 2 -c 2 (\)A)I-(c 2 (\)-c 2 (\))V®V, c?(A) > fc|(A) > 0. 

As will be more fully explained in Section 2.3, the null condition in this instance 

is! 

(N) Bi% n (\)UMiZmtn = 0, foralHeS 2 . 

With these essentials, we can now give a precise statement of the result 
to be shown in Section 3. 

Theorem 1.1. For every value of A £ M + such that the hyperbolicity 
condition (H) and the null condition (N) hold, the initial value problem for 
(1.9a)-(1.9c) with initial data 

du{0) £ Hi' 1 , k>9, 

has a unique global solution u £ H^(T) for every T > 0, provided that 

(1.12) £«-2(u(0)) exp [c(X)E 1 J 2 (u(0))] < e(A), 

ande(X) is sufficiently small, depending on A. The solution satisfies the bounds 

E K - 2 (u{t)) < e(A) and E K (u(t)) < 2E K (u(0))(tf^^ , 
for all t > 0. 

2. Preliminaries 

2.1. A word about invariants. In addition to their definitions as the 
elementary symmetric functions of the eigenvalues, the invariants 1(C) of 3 x 3 
matrix C are conveniently expressed as 

(2.1a) T 1 (C)=txC, J 2 (C) = i[(tr C) 2 - tr C 2 }, l 3 (C) = detC. 

The invariants of any two matrices which differ by a multiple of the identity 
are linearly related. Thus, if C' = z I + C, then 

(2.1b) 2i(C") = 3z + Ji(C) 

J 2 (C") = Sz^zX^+X^C) 

Z 3 (C") = z 3 + z 2 l 1 (C) + zl 2 (C)+l 3 (C), 

as can easily be verified by comparing eigenvalues. 



^Our earlier version was = 0, for all ^ £ S 2 and all fc. 
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In order to make sense of the coming assumptions about the stored energy 
function, it will be necessary to explore its dependence on the invariants of 
the strain matrix. When considering perturbations from an equilibrium, it is 
natural to introduce the invariants of the perturbation. We shall alternate 
between the strain matrix with which computations are most easily performed 
and its square root, the stretch matrix, in terms of which the results are most 
concisely expressed. The following table is meant to collect the notation to be 
used: 

(2.2a) l k (F T F) = ± k l k (F T F-X 2 I) = j fc k _ l23 

l k (VWF) = r k l k (VWF-XI) = s k 1 1 '"" > ' 

Thus, by (2.1b), i and r are linearly related to j and s, respectively. 

Through comparison of eigenvalues, the invariants i and r are seen to 
satisfy 

(2.2b) ii=r?-2r 2 , i 2 = r| - 2nr 3 , i 3 = r§. 

This transformation is invertible in a neighborhood of the equilibrium, by the 
implicit function theorem. From (2.2b), a short computation using (2.1b) 
produces 

(2.2c) j! = 2A Sl + s?-2s 2 , 

j 2 = 4A 2 s 2 + 2As!S 2 - 6As 3 + - 2s!S 3 , 

j 3 = 8A 3 s 3 + 4A 2 sis 3 + 2As 2 s 3 + s|, 
which is also invertible near the origin. 

Note. As a consequence of these formulae, it is possible to switch between 
any of the sets of variables i, j, r, s. In particular, by (2.1b), (2.2c), we can 
write 

(2.2d) a{F) = a{±) = t(A, s). 

2.2. The linear operator. In order that the linear operator L defined in 
(1.9a) be hyperbolic, we need to impose an ellipticity condition for the operator 
A in (1.9b). 

Lemma 2.1. Let A be the operator defined in (1.9b) from an isotropic 
stored energy function a. The following statements are equivalent: 

• A is elliptic. 

• The symbol of A, A 11 (£) = A l ^ m {\)^ m , is positive definite: 
(2.3a) <(A)^V>0, forallx^GS 2 . 

• The Legendre-Hadamard condition is satisfied: 

(2.3b) D 2 ct(\I + s x ® > 0, for all x,£ G S 2 . 

s=0 



856 



THOMAS C. SIDERIS 



• The positivity conditions 
(2.3c) th(A,0)>0 and n(A, 0) - At 2 (A, 0) > 

hold for the stored energy function r in (2.2d).l 

Proof. Condition (2.3a) is merely the definition of ellipticity. 

Thanks to (1.9b), the derivative in (2.3b) is equal to the expression in 
(2.3a), so the two conditions are equivalent. 

From (2.1a), it is a simple matter to check that when F = XI + e x £, 
x, £ £ S 2 , the invariants of F T F — X 2 I satisfy 

(2.4a) j 1 = 2Ae(x,£)+e 2 , j 2 = -A 2 e 2 |x A £| 2 , j 3 = 0. 

Combining (2.4a) and (2.2c), we find that S3 = and 

(2.4b) D eSi \ £=0 = (x, D £ s 2 \ e =o = 

£>?si| £=0 = ^x\xA £| 2 D 2 £ s 2 \ £=0 = -\\x A i\ 2 . 

So with (2. 2d), (2.4b), we can explicitly compute the derivative in (2.3b): 

(2.4c) D e V(A/ + ex®()|, =0 = D E 2 r(A 1 si,s 2 ,0)| £= o 

= r n (A,0) <x,£> 2 + 2i[ti(A,0)-At 2 (A,0)] |xA£| 2 , 

thereby showing the equivalence of (2.3c) and (2.3b). □ 

We shall now impose conditions (2.3a)-(2.3c). Considering (2.3c), we 
define positive constants ci(A) and c 2 (A), representing the propagation speeds, 
by 

(2.5a) c 2 (A) =r n (A,0), c§(A) = i[n(A, 0) - Ar 2 (A, 0)]. 

Following the standard linear theory, we shall assume that the reference con- 
figuration, A = 1, s = 0, is a stress- free state 

(2.5b) tl(1,0) = 0. 

It is necessary that the speeds remain distinct, and to be consistent with linear 
theory we assume that the bulk and shear moduli are positive: 

(2.5c) c 2 (A) - |cl(A) > and t|(A) > 0. 

From (1.10), (2.3a), (2.4c), and (2.5a) we see that the symbol matrix A{£) 
is positive definite and symmetric with the squared speeds as its eigenval- 
ues. The corresponding eigenspaces are the one-dimensional span of the unit 



4 Here and later on, subscripts for t indicate derivatives in S. 
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vector £ and its orthogonal complement. Consequently, we can write the sym- 
bol in spectral form 

(2.6a) A(£) = c?(A)£®£ + cl(A)[J-£®£], £ e S 2 , 

or in other words, 

(2.6b) A = c|(A)AJ + (c?(A) - c|(A))V (8) V. 

For each direction £ G S" 2 , we have two families of elementary plane wave 
solutions of Lu = 0, namely 

(2.7) Wi(£) = Kexpi/^O-dtA)*] :«,/?£!}, 

W 2 (£) = {??expi/3[(x,0-c 2 (A)t] : (rj,0 =0,/3gM}. 

These elementary solutions represent longitudinal and transverse waves prop- 
agating in the direction £ with speeds ci(A) and 02(A), respectively. 

2.3. TTie nonlinearity. Global existence requires a further nonresonance 
condition, the so-called null condition, linking the quadratic portion of the 
nonlinearity with the linear operator: the quadratic interaction of elementary 
waves of each wave family only produces waves in the other familyH This idea 
is expressed in the following: 

Definition 2.1. The quadratic nonlinearity N defined in (1.9c) is null with 
respect to the linear operator L defined in (1.9a) and (1.9b) provided that 

(2.8a) (u,N(v,w)) = 0, 

for all resonant triples 

(u, v, w) e W Q (£) x W Q (£) x W a (£), a = 1, 2. 

In terms of the coefficients BV m of the nonlinearity, we see by direct 
substitution that condition (2.8a) is equivalent to 

(2.8b) BZ k J\)^M^n = 0, for all (e5 2 , 

(2.8c) B^X^nfrSkiUU = 0, for all £, ^ G S 2 

with (£,?? (a) ) = 0. 

The following result relates this condition to the stored energy function 
and the stress tensor. In particular, it shows that condition (2.8c) for the 
transverse waves is redundant in the isotropic case. 

Lemma 2.2. Assume that the nonlinear quadratic form N defined in 
(1.9c) arises from an isotropic stored energy function a. The following state- 
ments are equivalent: 



5 A connection between plane waves and the null condition was first noted in [6]. 
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• The nonlinear form N is null with respect to L. 

• Condition (2.8b) holds. 

• The stored energy function satisfies 

(2.9a) D*a(\I + e £ ® Q\ £=0 = 0, forall£eS 2 . 

• The Piola-Kirchhoff stress tensor X(-F) = da(F)/dF satisfies 
(2.9b) tr£>f£(AJ + ef (g>OU=o£<S>£ = 0, for all £ G S 2 . 

• T/ie degeneracy condition 

(2.9c) r m (A,0)=0 
ZioWs /or r m (2. 2d). 

Proof. It is clear that conditions (2.8b), (2.9a), and (2.9b) are equivalent, 
while condition (2.8c) is equivalent to 

(2.10a) D £l D £2 D £:i a{\I + (e^ 1 ) + e 2 r/ 2 ) + e 3 ?? (3) ) ® €)| £a =o = 0, 
for all f.ijW G S 2 with (£,f/ a) > = 0. 

To see that this holds for all isotropic materials, take 

F = XI+ (eir/ 1 ) + e 2 r] {2) + e 3 ^ {3) ) ® £ 

with £,7/°) G S 2 and (£,7/^) = 0. It follows from (2.4a) with e = |e »7^| and 
x = e a rf a ^ je that now 

(2.10b) jx = le^W+e^W+eg^l 2 , 

j 2 = _A 2 |( £l?? ( 1 )+e 2?? ( 2 )+ £ 3r/ {3) ) A£| 2 , 
33 = 0. 

Note that (2.10b) is quadratic in £ a . Xherefore, regarding o~(i) as a function 
of j through (2.1b), we see that (2.10a) is true without further assumptions. 

Finally, to get (2.9c) let us consider (2.9a). Setting F = XI + e £ ® £, 
we get from (2.4a) that ji = 2Ae + e 2 , j'2 = J3 = 0, and so from (2.2c) it 
follows that si = e, S2 = S3 = 0. Switching variables, as in (2. 2d), we see that 
a(\I + e £®£) = r(A,e,0, 0), and the null condition (2.9a) for the longitudinal 
waves reduces to (2.9c). □ 

Note. When A = 1, the condition (2.9c) is complementary to John's gen- 
uine nonlinearity condition rm(l, 0) > which leads to formation of singular- 
ities in small spherically symmetric displacements [4]. 

Note. In terms of the Cauchy stress tensor, T(F) = det(F)~ 1 T,(F)F T , the 
null condition can be characterized as 

tr D 2 T(A/ + e £ ® OU=o £ <8> £ = 0, for all £ G S 2 . 
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2.4. Deciphering the conditions. This section will illustrate that our con- 
ditions can be satisfied for all values of A with physically realistic choices of 
the stored energy function. 

An expansion of the stored energy function can be based on the distor- 
tional strain matrix C = V F T F — |ri/, where, with the notation in (2.2a), 
n = tr VF T F is the dilational strain. (The distortional strain is essentially 
the linearization of the shear strain considered in [16].) Notice that tr C = 
and that C vanishes when F is a multiple of the identity. Since C can also be 
written as (V F T F — XI) — ^S]/, the remaining two invariants of C, z 2 = 2 2 (C) 
and Z3 = 23(C), can be expressed in terms of s, by (2.1b), and these are small 
for small displacements. Therefore, in terms of the variables, 

(2.11) zi = ri = A + isi, 

z 2 = 12(C) = s 2 - |sf, 

z 3 = 23(C) = s 3 - isis 2 + ^sf, 

the stored energy function has an expansion of the form 

(2.12) r(A, s) = /(zi) + g(z 1 )z 2 + M z i) z 3 + r(z 1 , z 2 , z 3 ), 

in which the remainder r(z 1; z 2 ,z 3 ) vanishes to first order in z 2 and z 3 : 

r(zi, 0, 0) = r 2 ( Zl , 0, 0) = r 3 (zi, 0, 0) = 0. 

For small s, it is clear that f(zi) is the dominant term in (2.12). 
The ellipticity conditions (2.5a) state that 

(2.13a) c\(X) = rn(A,0) = ^f'(X) - p(X), 

(2.13b) c|(A) = i-[ Tl (A,0)-Ar 2 (A,0)] = i-[I/'(A) - A 5 (A)]. 

Elimination of g from (2.13a) and (2.13b) leads to 

(2.14a) f"(X) - jf'(X) = 96(A), 

with 6(A) = c\ (A) — |c 2 (A), the bulk modulus which has been assumed to be 
positive in (2.5c). The ODE (2.14a) together with the initial condition (2.5b), 
which implies that /'(l) = 0, uniquely determines the function / (up to an 
inessential constant) in terms of the bulk modulus 

(2.14b) f(x) = 3£*-^b(y)dy. 

Having determined / in (2.14b) and therefore also g from (2.13b), the null 
condition (2.9c) uniquely determines the function h: 

= r m (A,0) = ±f"'(X) - Ig'(X) + ffc(A). 
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Starting from (2.12), the Cauchy stress for the homogeneous dilation 
<p(x) = Xx is T(AJ) = (/'(A)/3A 2 )/, in units of force per unit area. Thus, 
the pressure on the sphere \p\ = A is proportional to — /'(A)/A 2 . From (2.14a), 
positivity of the bulk modulus is seen to correspond to the monotonicity of 
this pressure in A. Experimental data for rubber-like materials suggest rapid 
growth for the pressure as A becomes small [11, p. 519]. It is clear from (2.14b), 
that as long as b(X) > 0(\~ p ) for A <C 1, with p > 2, the stored energy func- 
tion behaves correctly near the equilibrium: r(A, 0) — ► oo, as A — ► and as 
A — > oo. We emphasize that the null condition places no restriction on the 
bulk or shear moduli of the equilibrium. 

3. Proof of the existence theorem 

3.1. Commutation. When we recall the definitions of the coefficients in 
(1.9b), (1.9c), differentiation of the relations (1.7) yields 

(3.1a) 4 m (A) = A^(X)Q m Q Pj Q Si Q em , 

(3.1b) B £mnW = B Seri * WQaiQpjQ-ykQseQemQrin, 

for all proper orthogonal matrices Q. So A and B are isotropic tensors. 

Consider the one-parameter family of rotations generated by the Ug de- 
fined in (1.2b), 

Qe(s) = U e Q £ (s), Q e (0) = I. 
If tp : M. 3 — > M 3 , define the family of simultaneous rotations 
(3.2) T Q p(x) = Qp(Q T x). 

Under the hypotheses (1.7), this transformation leaves the equations of motion 
(1.8) invariant. 

From (1.9b), (3.1a) it follows that 

(3.3a) T Qe{s) [Au}=A[T Qe{s) u}, 

and in the same way, from (1.9c), (3.1b) that 

(3-3b) T Qe(s) N(u, v) = N(T Qe{s) u, T Qe(s) v), 

In particular, we conclude from (3.3a) and (3.3b) that (3.2) also leaves the 
truncated system (1.9a) invariant. 

The operators O defined in (1.2a) are generated by in the sense that 

(3.4a) n e u = D s T Qe ( s) u\ s=0 . 

Taking the derivatives of (3.3a) and (3.3b) at s = 0, we see that (3.4a) gives 
the commutation relations 



(3.4b) 



QiAu = AO,£U, 
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and 

(3.4c) n e N(u, v) = N(h e u, v) + N(u, fl t v). 

So from (3.4b) and (3.4c) any solution of the truncated linearized equation 
(1.9a) satisfies 

(3.4d) m t u = N(n e u, u) + N(u, n eu ). 

Next, define the one-parameter family of dilations 
(3.5) R s u(t,x) = s~ 1 u(st, sx). 

Since 

(3.6a) R s Lu = s- 2 L[R s u], 

and 

(3.6b) R s N(u, v) = s~ 2 N(R s u, R s v), 

the scaling (3.5) leaves the truncated equations (1.9a) invariant. 
The family (3.5) generates S defined in (1.3): 

(3.7a) Su(t,x) = D s R s u(t,x)\ s=1 . 

Hence, upon differentiation of (3.6a), (3.6b), we have that 

(3.7b) SLu = LSu - 2Lu, 

and 

(3.7c) SN(u, v) = N(Su, v) + N(u, Sv) - 2N(u, v). 

By (3.7b) and (3.7c) any solution of the linearized equations (1.9a) also satisfies 

(3.7d) LSu = N(Su,u) + N(u,Su). 

As a consequence of (3.4d) and (3.7d), we have 
PROPOSITION 3.1. For any solution u of (1.9a) in H^(T), 

(3.8) LT a u= N(T b u,T c u), 

b+c=a 

in which the sum extends over all ordered partitions of the sequence a, with 
\a\ < k — 1. 

Note. As a final remark in this section, we consider the commutation 
properties of the projections P a defined in (1.5a). Since Tq(x/r) = x/r, it 
follows that Qe(x/r) = 0. Thus, £le commutes with P a . Likewise d r (x/r) = 0, 
and so d r also commutes with P a . The projections do not commute with L, 
A, or V, however. 
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3.2. Pointwise estimates. The lemmas in this section use the decomposi- 
tion (1.1). Despite their simplicity, the next two form the heart of the decay 
estimates. 

Lemma 3.1. Let u e H\. For the linear operator A defined in (1.9b) and 
its corresponding propagation speeds c 2 a from (2.5a), 



P n 



< 



C 



Au(x) — c 2 a d 2 u(x 
for a = 1, 2. 

Proof. Thanks to (1.1), we may write 



\Vnu(x)\ + \Vu(x) 



jj = ^-5 r V — —d r (-^ Afl) v? - I 4 AO) 3 m u\ 
r z r \r z J m \r z 



from which it follows that 
(3.9) 



„2 r 



< 



c 



|VO«| + |Vu| 



by the commutation property (1.4). Recall (2.6a) which says that the symbol 
satisfies A(x/r) = X) a c^P a , so that the result is clear from (3.9). 

Lemma 3.2. Let u G fl?(T). T/ien for a = 1, 2, 



□ 



(3.10a) |cai-r||P Q Au(i,x)| < C 
and 

(3.10b) |c a i-r||P a <9 t Vu(i,z)| < C 



| VTu(t, x)\ + | Vu(t, x)\ + *|Lu(t, x)| 



| Vru(t, x)| + | Vu(t, x)|+i|Lu(t, x)\ 



Proof. The following easily verified identities appeared in [10]: 
(3.11a) 

(c 2 a t 2 - r 2 )Au(t, x) = c 2 a (td t - rd r )Su(t, x) 



— r 



Au(t, x) — (? a d 2 u{t, x 



c^t 2 Lu(t, x), 



(3.11b) 



(c a t — r)dtd r u(t, x) = (dt — c a d r )Su(t, x) + ^— ^-Au(x, t) 



+ 



Au(t, x) — c a d r u(t, x 



— tLu(t, x). 



The estimate (3.10a) is obtained from (3.11a) by dividing by c a t + r, 
applying P a , and using Lemma 3.1. 
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For the other inequality, we write 
(c a t - r)P a d t d e u(t,x) 



= (c a t - r) 



P a —d r d t u(t, x) + P a d t { d e - —d r ) u(t, x] 
r \ r 



The first term satisfies the desired estimate by (3.11b) and Lemma 3.1. The 
second term can be estimated after first writing 

(cat - r)d t u(x, t) = c a Su(t, x) - r[c a d r + d t ]u(x, t). 

x^ fx \ 

Then if we apply dt d r = — I —= A O ) to both sides, we find that the sec- 



ond term also has the right bound, since [ — A U] commutes with c a d r + dt- □ 

\r ) i 

Note. This lemma is the only spot where the scaling operator S is re- 
quired. A similar weighted estimate holds for d 2 u(t,x), but we have no need 
for it. 

The next result captures the manner in which the null condition will be 
useful in the course of the energy estimates. 

Proposition 3.2. Suppose that u,v,w £ H\. Assume that nonlinear 
form N satisfies the null condition (2.8a). LetM = {(a, (3, 7) 7^ (1, 1, 1), (2, 2, 2)} 
be the set of nonresonant indices. Then 

(3.12) \(u(x),N(v(x),w(x))}\ 

< —\u(x)\ [|VO a v(x)||Vw/(x)| + \Vn a w(x)\\Vv(x)\ 



\a\<l 



+ \V 2 v(x)\\Q a w(x)\ +\V 2 w(x)\\Q a v(x)\ 



+ Cj2\ p Mx)\ \Pf3V 2 v{x)\\P^w(x)\ + \P^ 2 w(x)\\P^v(x)\ 



Proof. Using the projections P a , we write B l e ^ n = £ B^P^ Pf P» k , 
so that 

(3.13) (u,N(v,w))= J2 BZ n PiP£ j P v a k u%{d m vid n w k ) 

a=l,2 

+ E B ZnP> 1 \p$ j d e d m vjpfd n W k + P» j d m vjP» k d e d n W k 

The second group of terms in (3.13) is estimated by the second expression 
on the right-hand side of (3.12). 



864 



THOMAS C. SIDERIS 



As in the proof of Lemma 3.1, the formula (1.1) enables the quantity 



u 



d e (d m vid n w k ) 



-d r {d r v j d r w k ) 



to be estimated by the first group of terms on the right-hand side of (3.12). 
The proof of (3.12) concludes since 

BZ n P2 l P^P V a kX -^=^ for all l ,j,k,a 

which follows by the null condition (2.8b) when a = 1 and by (2.8c) for 
a = 2. □ 

3.3. Sobolev inequalities. The following Sobolev-type inequalities with 
weights for the most part appeared implicitly in [10]. 

Lemma 3.3. For u G Co°(R 3 ) 3 , r = \x\, and p = \y\, 
(3.14a) r l / 2 \u{x)\ < C ^ l|VSl a u|| L 2, 

|a|<l 

r\u{x)\ < C ]T \\ d r^ a A%y\> r) ■ E H^Hl^i^)' 



(3.14b) 



|a|<l 



|o|<2 



(3.14c) r(c a t-r) l l 2 \u{x)\ < C £ \\{c a t - p)d r Sl a u\\ L * M > r) 

\a\<l 

+ C E ll^ an llL 2 (|j,|>r) ; 
|o|<2 

(3.14d) r(cat-r)|it(i)| < C £ ll(Ca* - p)d r n a u\\ LHly{ > r) 

\a\<l 

+ c ]T \\(c a t- p)n a u\\ L 2 ilyl > r) . 

\a\<2 

Proof. For the moment R(r) will denote any smooth and positive radial 
function. The coordinate on the unit sphere S 2 will be denoted by u and the 
surface measure will be dw. The proof of these inequalities begins with the 
bound 



(3.15) 



r 2+a R(r) 2+ P f \u{ruj)\ A dw < Cr 2+a f f 



R{p) 2+ P\d r u{pu;)\\u{pu)t 



+ R{p) l+ P\R'(p)\\u{p^ 



dpduj 
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< C 



\y\>r 



\R\ 2 \d r u\' A + \R'\ 2 \u {2 



1/2 



dy 



(7 \v\ 2a \R\ 

\J\y\>r 



2 "\n\^+fD\ u fdy) 



\ 1/2 

J 



At this point we pause to extract (3.14a) from (3.15) setting a = and 
R(r) = 1. With the use of the standard estimate \\u\\ L e < ||Vn|| L 2 in M 3 , we 
obtain 



(3.16) 



, \ 1/4 

|ti(rw)| 4 dcj) < C||Vu 



!l 2 - 



Making use of another basic Sobolev inequality 

(3.17) \u(x)\<CJ2 \\n a u(ru;)\\ Li{s2) , 

\a\<l 

we get (3.14a) from (3.16). 

We return to (3.15) in which we now take a = 2 for the remainder of the 
proof. We examine the last integral in (3.15), starting with the inequality 

IK r ^)llL6(S2) < C J2 ll^ au ( r ^)llL2 3 (52)ll u ( ra; )llL4 3 (52)- 
|o|<l 

From this we get 

(3.18) / \y\*\R\ 2 ( 1+ ®\u\*dy < Csup(Vl^(p)| 2+/3 / \u{puj)\ A doj 

J\y\>r p>r \ J S 2 

I \Rf\ti a u\ 2 dy. 

i i J \y\ >r 



\a\<l 



Putting (3.18) together with (3.15), we conclude that 
(3.19) 



r 4 \R(r 



)| 2+/3 J s Ju(ru;)\ 4 diu\ < C 



\>r 



\R\ 2 \d r u\ 2 + \R'\ 2 \u\ 2 

1/4 



1/4 



dy 



£ / \Rf\n a u\ 2 dy 
.i^-, J\v\>r 



Ja|<l 

The remaining inequalities are proved by combining (3.19) with (3.17). 
In order to get (3.14b), take R{r) = 1 again, while for (3.14c) and (3.14d) set 
R(r) = (c a t — r) (note that \R'\ < 1), and (3 = 0, 2 respectively. □ 
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The next result tunes these estimates to the demands of the proof by 
applying them to higher derivatives and also by removing the singularity near 
the origin. 

Proposition 3.3. Let u e H£(T), with X K (u(t)) < oo. 



(3.20a) {r^lT^x^KCE^iuit)), 
(3.20b) (r)\dT a u(t,x)\ < CE l J\u{t)), 
(3.20c) (r)(c a t - r) 1/2 \P a dT a u{t,x)\ 

<C E l J 2 {u(t))+X K (u{t)) 

(3.20d) {r){c a t-r)\P a dVT a u(t,x)\ <CX K {u{t)) 



\a\ + 2 < k, 
\a\ + 3 < k, 

\a\ + 3 < k, 
Id + 4 < k. 



Proof. Inequality (3.20a) and inequality (3.20b), for r > 1, follow by ap- 
plication of (3.14a) and (3.14b) to the indicated derivative and by the commu- 
tation property (1.4). On the other hand, inequality (3.20c) for r > 1, results 
from application of (3.14c) to P a dT a u and by the fact that d r and £1 politely 
commute with P a , as explained at the end of Section 3.1. Similarly, (3.20d), 
for r > 1, arises from (3.14d). (Since (3.14c) and (3.14d) delete the origin, we 
need not consider the lack of smoothness of P a u near the origin.) Thus, we 
need only prove (3.20a)-(3.20d) for r < 1. 

For r < 1, (3.20b) is an immediate consequence of the Sobolev embedding 

t3\ 



(3.21) 



H 2 (R 3 ) C 



To obtain the other inequalities, we define a cut-off function 
(3.22) CM 



1, ifr<l 
0, ifr>2. 



Consider (3.20a) for r < 1. By (3.21), 

\T a u(t,x)\ =((r)\r a u(t,x)\ < CJ2 l|V 6 (Cr a n)|| L2 

|6|<2 



< C[Ell 2 {u{t)) + \\T a u\\ LH ^ 2) ]. 



We finish by using (3.14a): 

ll ra ^llL 2 (|y|<2) < CWp^uU* 



\p- ll2 \\m\y\<2)<CEll 2 {u{t)). 



Next we turn to (3.20c), for r < 1. Note that there are constants C < C 
such that 



(3.23a) CCM EM " r > P " ^ CM<*>J < C"C(r/2) ^(c a t - r)P Q . 
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We again make use of (3.21) and (3.23a), 
(3.23b) (c a t - r) 1 ' 2 \P a dY a u(t, x)\ 

< C(t) 1/2 C(r)|9r a w(t,x)| 

< C(t) 1 / 2 J2 \\V b (0T a u)\\ L , 

\b\<2 

< c(i)V2 £ ||dv b r a u|| L2( |,|< 2) 

|6|<2 
P |fe|<2 

< CX K (u(t)) + CJ2 \\{cpt - p) 1/2 P(idT a u\\ L2(lyl < 2) . 

P 

To complete the proof of (3.20c), we must treat the last term above. By (3.14c), 
\\(c p t - P ) l,2 PpdT a u\\ L 2 { \ y \< 2) < \\ P (c p t - P ) l l 2 PpdT a u\\ L ^ Wp^Wv^k?) 

< C[E\l\u{t)) + X K {u{t))\. 

The proof of (3.20d) for r < 1 follows the same lines as (3.23b), except 
that the last step above is no longer necessary. □ 

Note. In the Lorentz invariant case, stronger versions of these inequalities 
hold which avoid the weighted norm X K altogether. 

3.4. Weighted I? -estimates. The goal of this section is to show that the 
weighted norm X K is controlled by the energy E l J 2 , for small solutions of (1.9a). 

Lemma 3.4. Let u e H^(T). Then 



(3.24) X 2 (u(t)) < C 



E l 2 /2 {u(t))+t\\Lu{t)\\ L , 



Proof. Recall definition (1.6): 

(3.25) X 2 {u(t)) = £ \\(c a t - r)P a dVu{t)\\ L 2. 

a 

In view of (3.10b), it is enough to consider the terms with two spatial deriva- 
tives. 

We introduce a smooth cut-off function along r = (c 2 t}/2, by defining 

(3.26) ^(t,r) = C(4r/<C2t)), 

using ( as defined in (3.22). On the support of if), the ratio (c a t — r) / (t) is 
uniformly bounded above and below, and so we have 

(3.27) CV(t, r)(t)I < 1>(t, r) £<c a t - r)P a < Cty(t, r)(t)I. 
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Thus by (3.27), the terms of interest have the estimate 
(3.28) £||<^t-r)P a V 2 u(t)|| La 

a 

< J2\\(cat-r}iPP a V 2 u(t)\\L^ 

a 

+ J2\\(c«t-r)(l-^P a V 2 u(t)\\Li 

a 

< (t)\\ipV 2 u(t)\\ L 2 + ]T || - r)(l - i>)P a V 2 u{t)\\ L 2. 

a 

First we will show that the integral over the support of ip is controlled by 
the elliptic operator A. From the explicit formula (2.6b), we have 

(3.29a) UAu{t)\\ 2 L2 = c\ J ^ 2 \Au{t)\ 2 dx 

+ 2c\(c\ - c\) J V 2 (An(t), V(V • u(t)))dx 

+ (c 2 -c 2 ) 2 J^ 2 \V{V-u(t))\ 2 dx. 

Since the first order spatial derivatives of ip are ©((i)" 1 ), integration by parts 
shows that 



(3.29b) J i(> 2 \Au(t)\ 2 dx 



Likewise, 
(3.29c) 



> ^||W m u(t)||| 2 -C(i)- 1 ||VVV*)llL 2 ||Vn(t)|| i2 . 



^ 2 (Au(t),V(V • u{t)))dx 



> ||^V(V • u{t))\\h ~ C{t)-^V 2 u{t)\\ L 4Vu{t)\\ L ,. 
From (3.29a), (3.29b), and (3.29c), we obtain, using (3.27) again, 
(3.29d) 

El l2 {u{t)) + {t)UAu{t)\\ L , 



J24(t)Ud e d m u(t)\\ L 2 < C 

£,m 

< c 



El /2 {u(t)) + ]T ||(c Q i - r)ipP a Au(t)\\ L 2 



On the support of 1 — ip, we note that the quantity (c a t — r)/r is bounded. 
So for the integrals over the support of 1 — ip, we use first (1.1) to pass from 
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an arbitrary second spatial derivative to radial derivatives, and then Lemma 
3.1 to pass to A, 



(3.30) 



(c a t - r)(l - il))P a did m u 



L 2 



< c 

< c 



El /2 (u{t)) + \\(c a t - r)(l - i,)P a d 2 r u{t)\\ L 2 



E l 2 /2 {u{t)) + c- 2 \\(c a t - r)(l - ^)P a Au{t)\\ L 2 
As a consequence of (3.25), (3.28), (3.29d), and (3.30), we have that 
X 2 {u(t)) < C7k 1/2 Kt))+^||M-r)P a ^)|| L2 

L a 

+ Y,W( c <*t-r)P a d t Vu(t)\\ L i 

a 

and so the lemma follows now from Lemma 3.2. 

Next we estimate the nonlinear terms on the right-hand side of (3.24). 
Lemma 3.5. Let u £ H^(T) be a solution o/(1.9a). Then 



□ 



with fx' 



Ey\u{t)) + X,,{u(t))E]J 2 (u(t)) + X^u(t))E 1 J 2 (u(t)) 



2 



+ 3. 



Proof. First applying Lemma 3.4 to T a u and summing over \a\ < fi — 2, 
we have that 

(3.31) X^u(t))<CE l J 2 {u(t)) + C Y, A\LT a u{t)\\ L 2. 

\a\<n-2 

Because of Proposition 3.1, we need only to estimate terms of the form 

t||v 2 r fe uvrs|| L 2, 

with \b+c\ = \a\ < fi—2. Plainly either |6| < m— 1 or |c| < m where m = [^-^-]- 
We note that the quantities {t)~ l {r)(c a t—r) are uniformly bounded below, 
so that 



(3.32) 



i^city^^Yi^t-^Pa. 



Thus we have the estimate 

t||v 2 r 6 nvr c u|| L 2 < cJ2 \\( r )( c at - r)p a v 2 r b uvr c u \\ L 2 



<c 



^ Q ||(r)(c Q t-r)P Q V 2 r b n|| L ^||Vr c U || L2 , \b\<m-l 
^ Q ||(c Q t-r)P Q V 2 r b n|| L2 ||(r)Vr c U || LOO , |c| < m. 
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In the first case, using (3.20d) we get the upper bound 

X m+ z{u{t))E l 'X{u{t)), 
and in the second case, by (3.20b) we get the upper bound 

^(«(t))^ 3 («(t)). □ 

The following bootstrap argument completes the estimation of X K . Notice 
that the smallness condition is imposed only for X K -2- 

Proposition 3.4. Let u e Hfi(T), k > 8, be a solution of (1.9a). // 
Entity)), n = k — 2, remains sufficiently small for < t < T, then 

(3.33a) < CE l J 2 (u{t)), 

(3.33b) X K (u(t)) < CE\l\u{t))[l + E]J\u{t)))- 

Proof. Since we have fj, > 6, it follows that = [^-} + 3 < fi. Thus, by 
Lemma 3.5, 



xMt)) < c 



and so we see that for E^ (u(t)) small enough, the bound (3.33a) holds. 

Since k > 8, it is true that k' = [^^] +3 < k — 2 = jjl. So again by Lemma 
3.5, we may write 

X K {u(t)) <C?[^/ 2 (u(t)) + A: M («(t))^/ 2 (u(t)) + A: K (u(t))Sj/ 2 (u(t)) . 

If E]J 2 {u{t)) is small, then this implies that 

A- K (u(t))<Ci?y 2 (u(t))[l + ^Kt))]. 
Thus, we obtain (3.33b) from this and (3.33a). □ 

3.5. Energy estimates. Suppose that u{t) G i?p(T) is a local solution of 
(1.9a). Let fi = n — 2, and assume that e(A), and hence E^ufo)), is small 
enough in (1.12) so that Proposition 3.4 holds at least for a short time. We 
suppose that T is the largest such time with this property. 

Proposition (3.1) and the energy method yield 

(3.34) K(u{t))= ]T ]T J(d t r a u,N(r b u,r c u)}dx. 

\a\<K~l b+c=a 



E]J\u{t)) + Ey\u{t))X,{u{t)) 
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Terms in (3.34) with b = a or c = a, |a| = k — 1, require special attention. 
Thanks to (1.10) they can be written as follows: 

(3.35) J(d t r a u,N(T a u,u)}dx = B^ n f d t r a U %(d m T a ^d n u k )dx 

= -B^ n Jd t d e T a u i a m T a vPa n u k dx 

2dt imn 



+ 2 



J d l T a u i d m T a u j d n u k dx 

d f r a u l d m r a u>d t d n u k dx. 



If we set 



N(u,v,w) = B^uVV, 



and 



(3.36a) E K (u(t)) = E K (u(t)) + ^ f N(T a u,T a u,u)dx, 

|0|<K-1 

then (1.11), (3.34) and (3.35) show that 

(3.36b) E' K (u(t)) Yl Yl f(dtT a u,N(T b u,T c u))dx 

|o|<K-l b + c=a 

+ J N(d t u,r a u ,r a u)dx. 



|o|<K-l' 

From (3.36b), we get 



E' K {u{t))<C J2 \\dT a u\\ L 2\\dVT b uVT c u\\ L 2. 



|o|<«— 1 b + c = a 



In the sum, we may assume that either \b\ < d — 1 or |c| < d, where d = [k/2]. 
As in the proof of Lemma 3.5, we get 



K«t)) < C(t)^El'\u{t)) 



X K {u(t))E l J^{u(t)), \c\ < d 
X d+ z{u(t))El /2 {u(t)), |6|<d-l. 



Since k > 9, we have that d + 3<fi = K — 2, and so by Proposition 3.4 we 
obtain 

(3.37) E' K {u(t)) < C(tr l E]/\u{t))E K {u(t)). 

(This estimate implies almost global existence, [5], [10].) 
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The next and decisive step is to show that the lower energy En(u(t)) 
remains small. Here finally the null condition enters the fray. From (3.36a), 
(3.36b), we have 



in which we again use the cut-off defined in (3.26). We shall focus on the terms 
involving the quadratic form N since the ones with the trilinear form N are 
similar. 

The integrals in (3.38a) over the support of ip where, it will be recalled 
the ratio (t)~ 1 (c a t — r) is bounded below, are estimated using (3.20c) and 
Proposition 3.4 



(3.38b) J ^(d t T a u,N(T b u,T c u))dx 

< c\\d t r a u\\ L 2\\^v 2 r b uvr c u \\ L 2 

< C{t)-^E]J\u{t)) 

x ]T \\(c a t - r)P a V 2 r^|| L2 ||( C/3 t - r) 1 / 2 P f) VT c u\\ L o 



since |6| + 2 < \i and |c|+3<^ + 1<k. Observe that the lack of a weighted 
estimate for the derivative Vr c u in L 2 forces us to estimate it in L°° , here and 
later on as well. 

To the second group of integrals in (3.38a), we apply Proposition 3.2. This 
leads to a bound of the form 





\a\<n-l L ' 




< C(t)- 3 / 2 ^ 2 (n(t))^| 6 | +2 ( U (t))[^3(n(t)) + ^ C | +3 (n(t)) 

< ^(t)- 3 / 2 ^/ 2 ^^))^^^))^^^)) 

< C{t)^ 2 E,(u{t))El l2 {u(t)\ 
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i|</*- 



(3.38c) C J2 J2 | ll^r a n|| L2 Hr-^l - V)Vr fe+1 n Vr c u|| L 2 

|o|<«-l b + c=a I 

+ ||r _1 (l- ip)V 2 T b u r c+1 n|| L 2 

+ S/ |-Paftr a «||-P/jV 2 r b u||P 7 vr c u|d 



Since r > (i) on the support of 1 — ip, the first group of terms in (3.38c) 
is bounded by 

(3.38d) c(t)-v 2 ]T £ ||5 t r^|| L2 

\a\<fj,— 1 fe+c^a 

b,cj^a 

|vr b+1 «|| L 2||rvr c -u|| L oo + ||v 2 r 6 «|| L 2||r 1 / 2 r c+1 u||Lo 

Since |6| + 2 < // and |c| + 3 < \i + 1 < k, we see by (3.20a) and (3.20b) that 
this expression is also bounded by 

c7<r 3/2 wo)^ /2 K*))- 

Now we consider the nonresonant terms in (3.38c), which will be estimated 
using (3.20c) and Proposition 3.4. Note that 

l<C{t)-V 2 {r){cpt-r){ Cl t-r) l /\ if (3 + 7 . 

Thus, if (3 7^ 7, we have 

(3.38e) 

y \p a d t r a u\\PfsV 2 r b u\\p^vr c u \dx 

< C(t)-* J \P a d t r a u\(cpt-r}\P^V 2 r b u\(r)(c y t-r)i\P^Vr c u\dx 

< C{t)-l\\d t T a U \\ L 2\\{cpt - r)P f) V 2 T b u\\ L 2\\(r)(G r t - r)^P 7 Vr%|| L oc 

< C(t)-V 2 E,(u(t))Ell 2 {u(t)), 

since \a\ + 1 < //, |6| + 2 < /i, and |c| + 3 < n. 

In the remaining cases where a ^ (3, the weight (r)(c a t — r) l l 2 can be used 
in front of d t T a u and estimated in L°° by E l J 2 (u(t)), with the help of (3.20c). 

Putting everything together (3.38a)-(3.38e), we have shown that 

(3.39) < C7(t)- 3 / 2 ^(n(t))Py 2 ( U (t)). 

The smallness of the energy implies, in turn, the smallness of |Vu| by 
(3.20b). Thus, by possible further restriction of e(A), the modified energy in 
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(3.36a) is equivalent to the standard one, for small solutions, and from (3.37), 
(3.39) we arrive at the coupled pair of differential inequalities 

K( U (t)) < c{t)-'Eyy{t))E K { U {t% 

%{u{t)) < C{t)-^E,{u{t))El/\u{t)). 

Now, the global bounds given in Theorem 1.1 are immediate, provided e(A) is 
small enough. 
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